m 

o 

(N 



Black branes dual to striped phases 



Benjamin Withers 



Q Centre for Particle Theory and Department of Mathematical Sciences 

^^ Durham University, South Road, Durham, DHl 3LE, U.K. 






> 



Abstract 



(^ We construct inhomogeneous charged black branes in AdS, holograph- 

^ ically dual to a phase at finite chemical potential with spontaneously 

j> broken translation invariance in one direction. These are obtained nu- 

S^ merically, solving PDEs for the fully backreacted system. Fixing the 

j^ periodicity scale, we find a second order phase transition to the inhomo- 

geneous phase. We comment on the properties of the state emerging at 
low temperatures. For some models we demonstrate the existence of a 
branch of striped solutions but no continuous phase transition. 



1 Introduction 

An approach to applied holography is to pick a gravitational bulk theory, either in 
a consistent truncation of supergravity or a phenomenological model, and search for 
solutions which describe various equilibrium phases of the dual field theory. Guided 
by condensed matter phenomena this strategy has proven to be a fruitful approach 
and has driven the construction of a variety of novel black brane solutions. 

One well studied class of solutions are those dual to holographic superfluids. In the 
simplest case a charged scalar field is added to the bulk model resulting in a branch of 
black branes with scalar hair, emerging from the charged Reissner-Nordstrom (RN) 
solution [1, 2, 3]. If this branch is thermodynamically preferred then the dual field 
theory will spontaneously break its global U(l) symmetry. The nature of the phase 
transition to the broken phase, if indeed there is one at all, depends strongly on 
the details of the bulk model away from the threshold temperature [4] . Such phases 
are not restricted to phenomenological models, existing in consistent truncations of 
supergravity Kaluza-Klein reductions [5, 6]. 

Another class of phases are those in which translational invariance of an electri- 
cally charged phase is spontaneously broken [7, 8, 9, 10]. Such phases are also seen 
for magnetically charged black branes [11, 12] and occur for a wide range of bulk 
models, recent examples include massive gravity [13] and Einstein-Maxwell-dilaton 
systems [14].^ In other examples, the instability can result in homogeneous but non- 
translationally invariant stationary black brane solutions, which can be constructed 
by solving ODEs [19, 20, 21]. In the case of [20, 21], the phase transition was found 
to be second order, leading to a modulated zero entropy state at zero temperature. 

In this paper we will focus on models of the type analysed linearly in [9], which 
break translations resulting in a non-homogeneous phase. The key ingredient in this 
model is a parity violating coupling with a pseudo-scalar field. A zero mode analysis 
indicates that this term is responsible for a new phase emanating from an electrically 
charged black brane, displaying both charge and current density waves. Moving 
beyond the linear analysis in these cases is desirable, principally to determine if there 
is a phase transition. Once a solution has been obtained one can begin to address 
issues such as the nature of the ground state, the number of solutions which can exist 
at a given temperature, the thermodynamically preferred spatial scale, transport 
properties, etc. 

In this spirit, we set about the construction of cohomogeneity-two charged black 



^See also [15, 16, 17, 18] for spatial modulation in a probe brane setting. 



brane solutions with AdS4 asymptotics. Whilst the linear analysis performed in [9] 
has the same cohomogeneity (hence 'stripes'), it does so without loss of generality; the 
zero modes may be simply superposed in different directions leading to solutions with 
lower symmetry. Thus we should expect cohomogeneity-three branches of solution 
emerging from the threshold temperature for linear stability, and ultimately these 
may turn out to dominate the ensemble. Here, for technical simplicity, we choose to 
focus on the special case of a cohomogeneity-two family of solutions, translationally 
invariant in one boundary direction. This requires the solution of a 2d PDE boundary 
value problem involving the Einstein equations. Our approach will utilise elegant 
methods [22, 23, 24] for rendering a problem such as this elliptic. In a technically 
similar but physically distinct scenario these methods have been employed in the 
construction of a variety of holographic lattices forced by UV boundary conditions [25, 
26, 27]. 

We note that solutions of the type constructed here have previously been sought 
in the work [28] employing a different numerical method and finding a first-order 
phase transition at temperatures lower than the threshold of linear instability. 

The organisation and key results of this paper are as follows. In sections 2 and 
3 we specify the bulk model, ansatz and numerical method. In section 4 we present 
a single striped solution. In section 5 we examine the temperature dependence at 
fixed periodicity, 27Cfi/k, finding second order phase transitions and evidence of a 
zero entropy inhomogeneous state at zero temperature. In section 6 we present an 
example of a different model, identical in linear theory, which does not exhibit a con- 
tinuous phase transition. That is, as with the holographic superfluid constructions, 
the behaviour away from Tc is strongly dependent on the details of the bulk model. 
We conclude in section 7. 

Note added: While this paper was being finalised [29] appeared, also finding sec- 
ond order phase transitions. 

2 The bulk theory 

We consider the bulk action 

^b.ik = /rfVv/=^(i?-^(90f-^F^-n0))-/^i^AF, (2.1) 



where F = dA and we have set IG-ttG = 1.^ This action will be supplemented later by 
Gibbons-Hawking and boundary counter-terms. We choose to work with functions 
r(0), V^(0) and '(9(0) which admit electrically charged RN black brane solutions with 
an asymptotic {L = 1) AdS4, when = 0. With this condition we may conveniently 
parameterise the functions in the vicinity of the phase transition in a small ex- 
pansion T = f + O(0)^ V = V + 0(0)^ and ^ = 1^ + 0(0)^ where f = 1- f 
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V = —6 + |r?7,^0^ and -d = ^%0. Note that such a description will only hold in the 
vicinity of the critical temperature and will not suffice for example, in determining 
the order of the phase transitions. 

Throughout this paper we will consider only the case m^ = —2, where we take 
the operator dual to the scalar 0, 0<^, to have dimension-2, setting its dimension-1 
source to zero. For concreteness we will consider two phenomenological models, 

r = sech(v^0), V = -6cosh(0/y3), ^ = ^ tanh(y30). (2.2) 

6v2 

r = f, V = V, ^ = ^. (2.3) 

The specific case of (2.3) with n = and ci = 4.5 was considered in [28]. The non- 
linear model (2.2) when ci = 6a/2 appears in a consistent truncation of a reduction 
on an arbitrary Sasaki-Einstein space, SEj [30, 31, 32]. Here we consider a phe- 
nomenological modification of this model by dialling the value of Ci away from 6-\/2- 
In practice we increase ci a small amount to ci = 9.9. This has the advantage of 
raising the critical temperature for improved numerical accuracy. 

3 Ansatz and implementation 

To construct the striped black branes we retain diffeomorphism invariance in the x 
and z directions, adopting the following stationary ansatz, 

ds^ = —i-Tf{z)dt^ + Z-—+X{dx + -fdzY + Y{dy + (3dtY) (3.1) 
^ V /(^) / 

A = Adt + Bdy (3.2) 

where T, Z, X, Y, 7, /3, A, B and are functions both of the AdS radial coordinate, 
z, and a single spatial boundary direction x. The ansatz is translationally invariant 
in the remaining spatial boundary coordinate, y. The function f{z) = {1 — z){l + 



^Note that to make contact with the hnearised instabihties of [9], we can identify T//ihcrG = 
■^T/fithcic as well as fc//ihcrc — ^fc/Mthcro- We have kept an equivalent linearised parameterisation 

of t(0), V{(I)) and -dicj)) such that n, Cihere = "-, Cuhere- 



z + z^ — fi'^z^/A) conveniently factors out the electrically charged RN solution at 
Z = T = X = Y = 1 and -f = P = B = (j) = and A = ij{l - z). This RN solution 
will constitute the normal phase of our system, with temperature T = {12 — fi"^) / (IGtc) 
and entropy density s = Air. We have chosen coordinates such that z = 1 gives the 
position of the outer horizon. 

3.1 A regular radial coordinate 

For presentational purposes the Schwarzschild-like 2;-coordinate is efficient, and we 
find that it is also sufficient for numerical purposes, however, by switching to a radial 
coordinate regular on the horizon we achieve improved numerical accuracy at low 
temperatures. Consequently in what follows we will adopt the radial coordinate r, 
defined by, 

{l-rf = l-z. (3.3) 

The black brane horizon is at r = z = 1 and the boundary at r = z = 0. 

We seek solutions with a regular horizon. This requires that the functions Z, T,X,Y, 
(],A,B, and are smooth functions of (1 — r)^ there, whilst the dxdr cross term, 
7 leads with (1 — r), but such that 7/(1 — r) is a smooth function of (1 — r)^. In 
addition, we choose the Killing vector dt to be null on the horizon; accordingly both 
A and /3 are taken to vanish there. '^ These conditions are imposed by factoring out 
appropriate powers of (1 — r) and using Neumann boundary conditions, for example, 
drB = at r = 1. Finally we have the additional Dirichlet condition T = Z at the 
horizon, rendering the surface gravity constant, determined by /(I). In particular 
the temperature of the inhomogeneous solutions will be the same as RN at the same 
/i, i.e. T = (12-;u2)/(167r). 

At the boundary we fix the boundary metric to be Minkowski and turn off all 
sources. Specifically we impose the Dirichlet boundary conditions Z = T = X = 
Y = l, P = 'y = B = (p = where we have now defined (p = np and A = fi at r = 0. 

These boundary conditions lead to a two parameter family of solutions charac- 
terised by the dimensionless temperature T/fi and periodicity k/fi. In particular 
there is no other freely specifiable data and consequently any spatial modulation will 
be spontaneous. 



•^We could have also allowed for (spatially modulated or homogeneous) boosted black brane 
solutions by setting /3 to a non- vanishing constant on the horizon, together with a suitable modified 
regularity condition for A. This would correspond to a Killing horizon for a linear combination of 
the Killing vectors dt and dy, and extend the family of solutions by one additional parameter. 



3.2 Implementation 

We may formulate the equations of motion resulting from (2.1) as an elliptic, bound- 
ary value problem by utilising the methods [22, 23, 24]. In short, the problem may be 
rendered elliptic by considering the Harmonic Einstein equations, obtained by taking 
the field equations of (2.1) and making the replacement, 

Rab — ^ Rab = Rab — "^ (a^b) (3.4) 

introducing the vector field ^" = g"^'^(X'^^ — r,'^^), where T is the Levi-Civita con- 
nection of g and F is the Levi-Civita connection of a reference metric g. Here, for 
simplicity, we take g to be the metric of the RN solution at the same value of /i. To 
ensure that we have a solution to the equations of motion of (2.1) we must verify 
that the vector ^ is zero to the required numerical accuracy. 

The resulting nine coupled Harmonic Einstein, Maxwell and scalar field equations 
are solved numerically using a Newton- Raphson method. The fields and their deriva- 
tives are represented using a spectral collocation method.^ We employ Chebyshev 
polynomial interpolants on an irregular grid for the r coordinate containing A^ + 1 
points in the interval [0, 1], and Fourier interpolants on a regular grid for the x coor- 
dinate containing A^ points in the interval [— vr/fc, ir/k) with x = —n/k and x = n/k 
identified. It will be convenient to denote quantities averaged over the period with a 

bar, i.e. 

k n^^ 

F{x)dx. (3.5) 



27r j^^/k 

We verify that our solutions are not Ricci-Solitons {i.e. that they are solutions to 
the equations of motion of (2.1)) by monitoring the field |.^^|, ensuring |^^| < 10~^° 
is satisfied at all grid points for the solutions presented in this paper. In practice 
higher values of A^ are required to achieve this condition at low temperatures, and 
so for the majority of solutions presented in this paper we achieve values which are 
smaller than this limit by several orders of magnitude. Tests of the numerics are 
presented in appendix A, in particular, in section A.l we verify that |^^| converges 
to zero exponentially with A^. 

As a further check, we have also constructed some of these solutions using fourth- 
order finite differencing, finding consistent results. 



*See [33] for a pedagogical exposition of these methods. 



4 A representative solution 

In this section we present a single representative striped black brane solution, for the 
model (2.2) with ci = 9.9. As an initial check of our numerics, we verify that the in- 
homogeneous phase meets the RN branch along a 'bell curve' of critical temperatures 
with good agreement with the n = 36, Ci ~ 9.9 linear instability case constructed in 
[9]. 

We focus here on a single solution in the one-parameter family of solutions emerg- 
ing from the linear bell-curve at k/fi = l.l/v2, for which the critical temperature is 
given by {T/fi)c — 0.0236. We select a solution at a temperature of T/ fi = 0.0205 ^ 
0.87(T//i)c. It is thermodynamically preferred over the RN solution. Bulk plots of 
the scalar and metric function Z for this solution are presented in figure 1. 

For this solution the Ricci scalar is shown in the left panel of figure 2, changing 
from i? = — 12 in the asymptotic AdS4 and becoming spatially modulated on the 
horizon. For comparison in the right panel of figure 2 we show the Ricci scalar for a 
solution on the same branch at much lower temperatures, T/Tc ~ 0.04, exhibiting a 
much larger variation on the horizon. 




Figure 1: Bulk profiles for the pseudo-scalar and a metric function, Z, for represen- 
tative striped black brane solution in the model (2.2), with ci = 9.9 and k/fx = 1.1/a/2 
at T/yU ^ 0.87(T/yu)c. The r-coordinate runs from r = on the boundary to r = 1 
at the horizon. 




Figure 2: Ricci scalar in the bulk. Left: The solution shown in figure 1, at at 
T//i ^ 0.87(T//x)c. Right: The same model and branch but at T//i ^ 0.04(T//x)c. 



5 A branch at fixed /c//i 

In this section we fix the periodicity of the solution, k/fi and examine the free energy 
as a function of temperature. As in section 4 we employ the model (2.2), with ci = 9.9 
and we focus on the one-parameter family at k/jj = l.l/-\/2. 

We find a second order phase transition at T/fi = {T/fi)c. The free energy for 
this branch, and its difference with that of RN is shown in the panels of figure 3. 
We find that the spatially modulated branch is thermodynamically preferred at all 
temperatures where it exists. To show that this corresponds to a second order phase 
transition, in figure 4 we plot the entropy density in the vicinity of {T/fi)c. According 
to the first law (which we test in section A. 2) the entropy s = — ||J at fixed /i, k. Hence 
we see that the kink in the plot of s/fi'^ at the critical temperature indicates a second 
order phase transition. 

An expression for the energy density is presented in appendix B. 



5.1 Lower temperatures 

For the data presented in figure 3, which we have obtained down to T/Tc = 0.1, we 
have ensured that the free energy has converged with A^ to sufficient precision, and 
that the first law test is satisfied (see section A. 2). However, the main source of error 
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Figure 3: Free energy for a fixed k/fi = 1.1/v^ branch of the model (2.2), with 
Ci = 9.9, indicating a continuous phase transition from the RN branch (black). 
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Figure 4: Entropy density for a fixed k/fj. = l.l/-\/2 branch of the model (2.2), with 
Ci = 9.9, in the vicinity of the phase transition from RN (black). Via the first law, the 
kink in this plot illustrates that the continuous phase transition at (T/fj,)c is second 
order. 



in this procedure is obtaining the energy density.^ If we are not interested in the 
free energy, then, we can go to slightly lower temperatures without compromising 
accuracy in any of the other observables, for example, in the entropy density s. 

The entropy density and the charge density are presented in the upper panels of 
figure 5 down to T/Tc ~ 0.03. In particular we see that the system appears to ap- 
proach an electrically charged, spatially modulated infrared geometry with vanishing 
entropy. The bottom two panels of figure 5 show the quantities v and j, defined as 
t>^ = ^ J (l)%)dx and j"^ = ^ J B?^sdx where 0(2) is the coefficient of z^ in a bound- 
ary expansion of and gives the scalar VEV, whilst -B(i) is the coefficient of 2; in a 



^This is presumably because it appears suppressed by z^ relative to the leading boundary be- 
haviour. 



boundary expansion of B and gives the current. 
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Figure 5: Various averaged quantities for a fixed k/jj. = l.l/-\/2 branch of the model 
(2.2), with Ci = 9.9, indicating an approach to an electrically charged, spatially 
modulated zero-temperature state with vanishing entropy. The quantities v and j 
are defined in the text. The RN solution is in black. 



6 The model (2.3) 

In this section we comment on a different bulk model (2.3), where we find that a 
continuous phase transition (at fixed k/fi) is absent at our parameter choice. For 
comparison with the results of section 5 we work with ci = 9.9, n = 36 and study 
the fixed k/ n = 1.1/a/2. This model is obtained simply by expanding the functions 
V, r, i} of the model (2.2) to first non-trivial order in (p. Consequently, a linear analysis 
yields identical results, e.g. for the critical temperature {T / n)^ — 0.0236. 

Perhaps unsurprisingly, we find drastically different behaviour to the model (2.2) 
once we move away from {T/fi)c. In particular we find that the branch of solutions 
emerging from the RN branch exist initially at T/fi > {T/fi)c and are subdominant. 
This is illustrated by plots of the free energy in figure 6, near the phase transition. 
For completeness we show the entropy in figure 7. 
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Figure 6: Free energy for the fixed k/fj. = l.l/-\/2 branch in the model (2.3) at 
Ci = 9.9, n = 36, demonstrating the absence of a continuous phase transition from 
the RN branch (black). 




Figure 7: Entropy for the fixed fc//i = l.l/v2 branch in the model (2.3) at ci = 
9.9, n = 36. The RN entropy is shown in black. 

7 Discussion 

In this paper we have constructed new inhomogeneous black brane solutions in AdS, 
describing spatially modulated phases of a dual field theory at finite chemical poten- 
tial. These modulated phases arise in the dual CFTs as new branches that sponta- 
neously break translation invariance. We studied branches of striped black branes at 
fixed k/fi in two models finding a range of possible behaviours. 

In the first model (section 5), as the system is cooled we find a second order 
phase transition to the striped phase. This is our main result. Additionally, at low 
temperatures we have found indications that a zero entropy, electrically charged, 
inhomogeneous state is being approached. We note the similar behaviour of these 
striped solutions to the helical case [20, 21] exhibiting a zero entropy, homogeneous 
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modulated state at T = 0. It would be interesting to go to lower temperatures and 
investigate this 'ground state' more thoroughly, as well as to construct it directly at 
zero temperature. 

In the second model (section 6), taken to be identical to the first in linear theory, 
we illustrated the existence of a branch of striped solutions at higher free energy 
than the RN branch. Moreover, these solutions exist at higher temperatures than the 
point at which the two branches meet, {T/fi)^. There is therefore no continuous phase 
transition. The natural question is whether a first order transition occurs from RN 
for T/yU > {T/fi)c e.g. with a swallow-tail structure. We note that this possibility 
is distinct from the first order transitions presented in [28], where the transition 
temperature was found to be lower than (T/yu)^. If no first order transition occurs 
then the situation becomes reminiscent of the 'retrograde' branches of holographic 
superconductors[34, 35, 4]. Irrespective of what happens for this particular model 
further along the branch, we have shown that a continuous phase transition need not 
occur and that, as expected, the nonlinear details of the model are important away 
from the critical temperature. 

Finally, it would be interesting to relax the two main assumptions adopted in this 
work, namely, homogeneity in the second spatial direction, y, and the restriction of 
fixed k/fi. In particular, it may be that solutions with lower symmetry dominate the 
ensemble, and that the thermodynamically preferred periodicity changes with the 
temperature, as in [20, 21]. 
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A Numerical tests 

A.l Testing convergence with A^ 

Using the solution presented in section 4 we will illustrate the convergence properties 
as we scale the number of points in the grid, A^(A^ + !)• As we are using spectral 
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methods we should find an exponential rate of convergence with A^. This is illustrated 
for the average of the charge density p in the left panel of figure 8. Crucially, we must 
ensure that the maximum value of |^^| also converges appropriately towards zero so 
that we have a solution to the Einstein equations - this convergence is indeed seen 
and is demonstrated in the right panel of figure 8. 
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Figure 8: Left: convergence of the charge density p with the grid size A^. Right: 
convergence of the maximum value of |^^| on the grid towards zero with increasing 
grid size A^. Dashed lines are the best linear fit. 



A. 2 Tests of the first law 

Combining the local free energy w{x) = e{x) — Ts{x) — fip{x) with the first law 
de{x) = Tds{x) + fidp{x) we conclude that along our fixed k/fi branch of solutions the 
quantity W = -^pf- + s/p"^ must vanish. In particular, approximating the derivatives 
by the difference of two nearby solutions with parametric separation X = 6 T/fi, we 
should expect to find numerically that W = -|- 0(A). Indeed, in figure 9 we see that 
this condition is met along the set of solutions of section 5, where A ~ 10~^. 



B Free energy 

To calculate the free energy we follow [36, 32]. Including boundary terms, the action 
for this system is given by. 



S = S, 



bulk 



d^xV^{K - 4 



(B.l) 



where Shuik is given by (2.1), hat is the induced metric on the boundary and Kab 
the extrinsic curvature. From this we may extract the renormalised boundary stress 



12 



-2.0 








' 


-2.5 






: 


-3.0 




, 


. . . • : 


\og,oW-3.5 


. . . 


■ 


. , 


-4.0 






': 


-4.5 


". 


, 


: 


-5.0 


- 




: 


-5.5 








0.005 


0.010 0.015 


0.020 






T//X 





Figure 9: Numerical check of the first law W for a set of solutions on a fixed k/^ 
branch. Note that because W is approximated here by differencing solutions with 
separation A = ST/fi, it is limited by A rather than N, and as required we find 



VT ~ A~ 10 



-3 



tensor T, 



ab 



SS 



^Ih^^ which is given in terms of the data appearing in the boundary 
expansion. From Tu we may identify the local energy density, 



ex 



2 



3T(3)(, 



(B.2) 



where T^'^\x) = |i9fT|^^p. This allows us to evaluate a local expression for the 
thermodynamic potential w{x) = e{x) —Ts{x) — fip{x). The entropy density is given 
in the usual way by s{x) = ^^ where a{x) is the area measure of the horizon. For 
the thermodynamical properties of the solutions presented in this paper we work with 
the averaged quantities, w,e,s,p, defined by (3.5). 
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